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Abstract 

Electromagnetic wave propagation through cold collision free plasma is studied using 
the nonlinear perturbation method. It is found that the equations can be reduced to 
the modified Kortweg-de Vries equation. 

1 Introduction 

An exciting and extremely active area of research investigation during the past years has 
been the study of solitons and the related issue of the construction of solutions to a wide 
class of nonlinear equations. The concept of solitons has now become umbiquitous in 
modern nonlinear science and indeed can be found in various branches of physics. In 
nonlinear wave propagation through continuous media, steepening of waves arises due 
to nonlinearities which is balanced by dissipative or dispersive effects. Exciting and im- 
portant discoveries were made in the nonlinear dynamics of dissipative and conservative 
systems. There are different methods to study nonlinear systems. The reductive pertur- 
bation method for the propagation of a slow modulation of a quasimonochromatic wave 
was first established by Taniuti and Washimi for the whistler wave in a cold plasma. This 
method was generalised to a wide class of nonlinear wave systems by Taniuty and Yajima. 
Kakutani and Ono [1], Kawahara and Taniuti [2] and Taniuti and Wei [3] have investigated 
the propagation of hydromagnetic waves through a cold collision free plasma using this 
reductive perturbation method. 

In the study of the asymptotic behaviour of nonlinear dispersive waves, Gardner and 
Morikava [4] were the first to introduce the scale transformation 

£ = £ a (x-vt), T = e (3 t. 

This scale transformation is called the Gardner-Morikawa [4] transformation. They com- 
bined this transformation with a perturbation expansion of the dependent variables so 
as to describe the nonliear asymptotic behaviour and in the process they arrived at the 
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Kortweg de-Vries [KdV] equation [5] which is a single tractable equation describing the 
asymptotic behaviour of a wave. This method has established a systematic way for the 
reduction of a fairly general nonlinear systems to a single tractable nonlinear equation de- 
scribing the far field behaviour. The reductive perturbation method was first established 
for the long wave approximation and then for the wave modulation problems. 

In the present work we study the propagation of electromagnetic waves through a cold 
collision free plasma by using a nonlinear reductive perturbation method. It is found 
that to the lowest order of perturbation the system of equations can be reduced to the 
modified Kortweg-de Vries equation (mKdV) [6]. In the case of steady state propagation 
this equation can be integrated to give a solution in terms of hyperbolic functions which 
exhibit solitary wave nature. 



2 Formulation of the problem 

When electromagnetic waves pass through a medium , the system gets perturbed. Since 
electrons are much lighter than ions, electrons respond much more rapidly to the fields 
and ion motion can be neglected. In the equation of momentum for cold plasma, no 
pressure term is present. Basic equations relevant to the present problem are the equations 
of motion of electron and the Maxwell's equations. Here we are interested only in the 
electronic motion. To obtain a single equation which incorporates weak nonlinear and 
weak dispersive effects, we employ the expansions of the dependent variables similar to 
that introduced by Nakata [7]. 

The equation of motion of an electron in an electromagnetic field is 



E + vx B 



dv e 
dt m 

Taking the leading order terms we get, 
dv e r 



(1) 



dt 



m 



E + vxB 



For convenience we take the displacement vector field S, which describes the direction and 
distance that the plasma has moved from the equilibrium. That is, 



_. dS ,_. . -=* 

— «+(-V)S. 



Therefore Eq.(|]) can be written as 

d 2 s 



dt 2 



e 
m 



E + 




(2) 



From Maxwell's equation 
and writing 



E = Eq exp i(kx — ut) 
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we obtain: 

kE = uB. 

Substituting Eq.(0) in Eq.pl) we can write: 



d 2 s 
dt 2 



E + 



k 



ds 

— - x E 
dt 



W, 



(3) 



(4) 



where W = — . Physically each electron is acted upon by an electric field that is parallel 

u> 

to the magnetic field so that there is no perpendicular component of motion that could 
be affected by the Lorentz force, that is v X B = [8]. Then Eq.(|l|) can be written as 

dv -> 
m— = —eE. 
dt 

then becomes 



Equation ( 

dt 2 

which can be put as: 
dS 



dv dS -A, 



dt 



v+(S x E)W. 



From Maxwell's equations, we have 
B = hqH, Vx£ 
dE 



dB 
~dt : 



V x H = J + e - 



dt 



Taking the time derivative of Eq.(^) we get, 

,d 2 E dv 



VxVx£ = 1/c 2 



dt* +fi0 d~f 



Equation (^) can then be written as 

2^2^ d 2 E d 2 S 
c 2 V 2 E = -^r + 



(5) 

(6) 
(7) 

(8) 
(9) 



dt 2 dt 2 ' 

Equations (a) and @ are systems of complicated nonlinear partial differential equations 
for E and S describing electromagnetic wave propagation through plasma. Let us seek 
a solution of these equations in the form of a Fourier expansion in harmonics of the 
fundamental E = e l ( fea:_a; *) as, 



E= J2 E n ' En i 



(10) 
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s= J2 s nEn - ( n ) 

ra=— oo 

Let us now consider one dimensional plane wave propagating along the x direction in 
the Cartesian coordinate system (x,y,z). All the physical quantities are assumed to be 
functions of one space coordinate x and time t. We now introduce the streching variables 
£ and r as, 

Z = e(x-Vt), 
t = e 3 t, 

where the velocity can be determined by the solvability condition of the above equations. 
E and S satisfy the following boundary conditions, 

E x ->0 except E° x = E cos 9, (12) 

£j-»0 except E° = E sm9, (13) 

E\ - 0, (14) 

as i->-oo, % = 0,1,2,3, ... (15) 

S x ^0 except S° = S cos 9, (16) 

Sj-fO except 5° = 5 o sin0, (17) 

SI - 0, (18) 

as £^ -oo, i = 0,1,2,3,... (19) 
The operators in terms of the streching variables can be written as 

dx 2 [ } ae 

For an appropriate choice of the coordinate system we can write E = (E x , E y ,0) and 
S = (S x ,S y ,0). Expressing the Fourier components of E and S in powers of a small 
parameter e 

oo oo 

S n = Y / e^S n J (x,t), E n = J2e j E?(x,t). 

j=0 j=0 
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Before proceeding to the nonlinear problem, it may be instructive to examine the 
dispersive relation in the linearized limit. Assuming a sinusoidal wave expi(kx — uit), 
where k, u are respectively the wave number and the frequency of the wave. Expanding 
the above coupled equations, 

d . 



dt 



moo 



S n = J2 S p xE q , 

p+q=n 



d 2 d 



d 2 , d 2l2 
+ 2ink-r n k 



dx 2 



dx 



Eq.@ gives the components of Si,™ as funtions of E\™, (s = x, y, z) 
The determinant of this system, A(n) is 



A(n) = inuj —n 2 ^ 2 u) 2 + n 2 s 2 + 7/x(l + a)sf 



where 



H = (1 + a<y)W, 7=1 



o 



UJ 



Eq 

5*0 



For n=l, A(l) is zero if oj satisfies the dispersion relation 



7 2 w 2 + fi 2 s 2 x + 7 /i(l + a)s 2 



0. 



From this we obtain V 
constants and that 



a 



(20) 



[1 + a 



-c. We assume that Sq = s and Eq = as are 



M n = Hq = for n + 0. 

The assumed conditions at infinity are, E^S™ — > for j = for all "n" except for 
(J, I n |) = (1, 1), where the limit is assumed to be a finite constant. For n = 1, A(l) = 
for j = 1. Under this condition the system has a nontrivial solution. But for n = 2, 3, 4, . . ., 
An ^ 0, we have the trivial solution. That is for j = 1 and n > 1, we get E 1 " = 5f = 0. 

For n = 0, A(0) = 0, we can choose E® = S® = 0. This completes the solution at order 
(l,n). 

For the next order, we can proceed in the same manner. The system will have a solution 
only if the determinant of the augmented matrix is zero . 
Now expanding the dependent variables as, 



S x — So + e S„ + e l St. + • • 



2 c2 



Sy Sy ~\- E Sy "j" £ Sy ~\~ ' ' 



S z = S® + £ l Sl+£ 2 S 2 + 



E x — E® + £^E\ + £ 2 E 2 + ■ 

Ey = Ey + £ Ey + ^Ey ~\~ " 
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E z = E° z +e 1 E 1 z +e 2 E 2 + ---. (21) 

Subtituting these expansions in Eqs.([|) and (|i~0|), then collecting and solving coefficients 
of different orders of e J for n = 1 with the boundary conditions given by Eqs.(|l^) to (17) 
we get: 
at order e° 

qO ipO qO pO _ n 
& x ti x - b x & z - U, 



da 



2 



0. 



^ 2 ^(^ + ^ ) = °; (22) 

at order e 1 



(l + a)5°^ = -y- ? ' 
Ocl 



(l + a)5^ = -y 



OH 



(El + 5i.) = 0, 

(El-aS l y ) = 0, (23) 

( J E 2 1 -a5 2 1 ) = 0; (24) 
at order e 2 

V m = S ' {El ~ aSl) ~ S * {E * ~ a ®> (25) 

= S° X (E 2 - aS 2 y ) - S° y (E 2 x - aSl) + (1 + a)^Sj, 

(E 2 x + S 2 x ) = 0, (26) 

gW-^) = _ 2y (i^£^ > (27) 

a£ c z or 
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8(E 2 -aS 2 



-IV 



(1 + a) 2 8S° Z 



d£ c 2 8t 



Solving for E 2 , S 2 and E 2 , S 2 , from Eqs.fl23|) and (^) we can get 



(E 2 y -aS 2 y ) = J -2V(l+af/c^dH, (E 2 -aS 2 z ) = J -2V(l+a) 2 / c 2 |^.(28) 



V 



dSl 



S° y {El-aS-i)-S»{E*-aSl). 



Of ei2 



Substituting for 

{El - aS 2 y ) and (E 2 - aS 2 ) 
from Eq.(p7D in Eq.fl25|) we get 



5 



-2V(l + a) 2 8 
c 2 8t 



S° z dC - s° z 



-2V(l + a) 2 8 c0 



S y J ^ S zdi ~ S° z J ^S° y dC - - — 



Or 
dSl 



« = V 



on ' 



(1 + a) 2 d£ 



Now introducing two new variables A and 9 defined by 
Sl = Acos6, S° z =Asm9, ^ = 5 sin^ i 



as £ — > oo. 



Equation (21) can be written as 

oi _ V gg 
(l + a)0f 

Now subtituting the value of and using the new variables Eq.(p6[) can be written as, 



cos 



8 f 8 C d 2 9 

9— / sin 6d!r — sin#— / cos9d!r = —ll—k. 
8t J or J 8£ 2 



Differentiating Eq.(28) with respect to £ and simplifying we obtain 
89 d 3 9~ 



8_ 



89 



8 2 9 89 
WW 



This can be integrated with respect to £ to give, 



89_ 8 3 9 
8r~ + ^W 



8 2 9 89 „ 
-To—dt,- 



8? 8^ 
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Fig. 1. Shows the variation of /(£) with respect to £. 



dO 

Multiplying throughout by — we get, 

OT 



30 d 3 1 fd0\ 2 80 



- 



dO 

Putting / = — , the above equation becomes 
df 3 „ 2 df d 3 f 

This equation is the modified Kortweg-de Vries (mKdV) equation. In the case of steady 
propagation of the wave this equation can be integrated to give a soliton solution 

/(C) = 2a sech (a () 

With £ = £ — At, A = constant, a 2 = — if and only if A > (fx > 0). 

A* 

Fig. 1 shows the variation of /(£) with respect to £. Figs. 2 and 3 show the variation 
of /(C) with respect to £ and r for different values of A. Since / = — , is defined as 

= arccos(l — 2sech 2 aC). 

It is seen that increases from to 2ir or decreases from to — 2tt according as a > or 

C 

a < o as C goes from — oo to oo, since is given by = J fd(. 



3 Conclusion 

Starting from the basic equations describing the propagation of electromagnetic waves 
through the cold plasma we have showed that the system of equations can be reduced to 
mKdV eqation. 
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Fig. 2. Show the variation of /(C) with respect to £ and r for A = 21. 



Fig. 3. Show the variation of /(£) with respect to £ and r for A = 30. 
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